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^ ■ Abstract 

The Mathematica toolkit AMBRE derives Mellin-Barnes (MB) representations for 
Feynman integrals in d = 4 — 2e dimensions. It may be applied for tadpoles as 
well as for multi-leg multi-loop scalar and tensor integrals. AMBRE uses a loop-by- 
loop approach and aims at lowest dimensions of the final MB representations. The 
present version of AMBRE works fine for planar Feynman diagrams. The output may 
be further processed by the package MB for the determination of its singularity 
structure in e. The AMBRE package contains various sample applications for Feynman 
integrals with up to six external particles and up to four loops. 



X 

5-( 



Contents 

1 Introduction 

2 Construction of Mellin-Barnes representations 

3 Using AMBRE Q 

4 One-loop integrals [ll| 

4.1 Example: the pentagon diagram of massive QED [ll 

4.2 Numerators QJ 

4.3 More masses [15 

4.4 More legs 

5 Multi-loop integrals: loop-by-loop integrations 

5.1 Example: two-loop planar box in massive QED [T7 

5.2 Special numerators 

5.3 Further examples: A three- loop planar box, a four-loop self-energy, and a 
two-loop pentagon Qj^ 

6 Tadpoles [21 

7 On-shell diagrams [22] 

8 Non-planar topologies 

9 Summary |24 
A AMBRE functions list [25 
References raj 



2 



List of Figures 



1 Massive QED pentagon diagram (examplel.nb, example2.nb) [12 

2 Massive QED one- loop box diagram (example3.nb) [14 

3 General one-loop vertex (example4.nb) [15 

4 Six-point scalar functions (example5.nb); left: massless case, right: massive 

case [16 

5 Massive two-loop planar QED box (section 5.1: example6.nb, section 5.2: 
example7.nb) [l_7 

6 The loop-by-loop iterative procedure (example8.nb) [20 

7 Massless topologies; left: four-loop two-point diagram (example9.nb), 

right: two-loop five-point diagram (examplelO.nb) [21 

8 Four- loop tadpoles with three massive lines; left: with a one- dimensional 
MB-representation (examplell.nb), right: with a six- dimensional 
MB-representation (examplel2.nb) [22 

9 On-shell self-energies; left: SE513m2 (examplel3.nb), right: SE311m 
(examplel4.nb) [23 

10 Non-planar massless vertex [23 



3 



1 Introduction 



Recently, Mellin-Barnes (MB) representations of Feynman integrals have been used ex- 
tensively in various phenomenological and theoretical studies of quantum field theory. In 
many applications, sometimes in quite sophisticated ones [1,2,3], the MB- integrals can 
be solved analytically. One also may merge knowledge of some analytical solutions given 
by MB-integrals with other methods, e.g. the differential equations approach, as demon- 
strated in [4]. An introduction to the subject with many examples may be found in the 
monographies [5,6]. A systematic derivation and numerical evaluation of MB-representa- 
tions for Feynman integrals with a (unpublished) Maple package was described in [7]. At 
the same time, the Mathematica program MB for the automatized analytic continuation 
of MB-integrals was published in [8]. With AMBRE, we deliver a Mathematica tool for 
the derivation of MB-integrals and their subsequent analytic continuation and numerical 
evaluation with MB. 

The article is organized as follows. In section 2 we introduce the formulae used for the MB- 
representation of a general Feynman integral. The basic features of AMBRE are described in 
section 3. One-loop examples are given in section 4. Section 5 describes the implementation 
of the loop-by-loop approach to multi-loop integrals. Examples with tadpoles and on-shell 
diagrams as well as problems related to non-planar topologies are discussed in sections 
6-8. A summary follows in section 9. In an appendix we list the Mathematica functions 
of AMBRE. 



2 Construction of Mellin-Barnes representations 



The backbone of the procedure to build up MB-representations with AMBRE is the relation 



— - = — / daA a B~ a T(-ar)T(u + a), (1) 

(A + B) u 2mT(v) J K ' v h y ' 

—ioo 

where the integration contour separates the poles of the T-functions. 

The object to be evaluated by AMBRE is an L-loop Feynman integral in d = 4 — 2e 
dimensions with N internal lines with momenta and masses m,, and E external legs 
with momenta p e : 



c lT(k)] = 1 r d%...d d k L T{k) 

11 Wl (in*/*) L J {ql - ml)* . . . {qf - m?)"i . . . {q 2 N - m^)*" ' K ' 



1 Often one uses the additional normalization e £lEL ; we leave this to the later evaluation with 
the package MB [8] . 
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The numerator T(k) is a tensor in the integration variables: 

T(*) = l, (3) 

The momenta of the denominator functions di may be expressed by external and loop 
momenta: 



(4) 



d i = q 2 i -m 2 i = [y^auh- PA - m 2 = i^aak - ^2 (3 ie p e \ - 

\l=l ) \l=\ e=l / 

In the package AMBRE, in a first step the momentum integrals are replaced by Feynman 
parameter integrals: 

G L [T(k)] = n / r( ^ J j Rdxj x^Kl -X>) { F [ x) N v - dm P L(T) 



o i =1 

(5) 



with 



N 



N v = Y,"i- (6) 



i=i 



The two functions C/ and F are characteristics of the topology of the Feynman integral. 
One may derive them from 



N 

M = Y / x i (q 2 -m 2 ) = kMk-2kQ + J, (7) 
i=i 

where M w = Y^?=i ®U'®uXi, and Qi = E^Ii unPi^i, and J = E J J Ii(-^ 2 - w-)^; namely: 



£/(x)=det(M), (8) 
F(x) = -det(M) J + QMQ. (9) 

The £7 and F as well as M = det(M) M~ l are polynomials in x, and so are the numerator 
functions Pl(T) in (2) for scalar and vector integrals: 

Pl(1) = 1, (10) 
Pl(A?) = E Af«'Q?- (11) 
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Tensors of higher degree depend additionally on the diagonalizing rotation V for Af, 



N diag = {a u ...,a L ) = (V'YMV- 1 , 



(12) 



and become non-polynomial in x. As an example, we quote here the case of an L-loop 
integral with a tensor of degree two: 



i=i 



[MiiQi] a [MviQif - 



r ( N ^-i L - 1 ) UF (v l -Y(v^) g 



a/3 



r (n v - |l) 

The formulae simplify considerably for one-loop integrals: 



on 



(13) 



N 



U = M = M = det(M) — V — J2 x i = !. 
F = -UJ + Q 2 = J^PiPj ~ Pi + mfaiXj = ^fijXiXj. 
Then, the tensor factors P\{T) in (5) will become: 



(14) 
(15) 



^i(l) = l, 

Af 



N 

E 

i=l j=l 



i=l 
N 



p ( AT d I 



r (^-|) 



a/3 

F ^— , etc., 



(16) 
(17) 

(18) 



with being the so-called chords introduced in (4). For the general case Pi(T) see 
section 4.2. 

One now has to perform the x-integrations. In AMBRE, we will do this by the following 
simple formula: 



1 Af / 
n i=l V j 



Xi = 



r(qi + --- + q N )' 



(19) 



From the above text it is evident that the integrand of (5) contains besides simple sums 
of monomials Jlj X T a ls° different structures. This is due to the appearance of the factors 
U(x) and F(x). Beginning with two-loop tensor integrals, one faces additionally a com- 
plicated dependence of P(T) on x for higher rank tensors T due to the appearence of V 
and a, see (13). 
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For this reason, the present version of AMBRE is restricted to scalar and vector integrals 
and/or to one-loop integrals. In these cases one may rewrite F(x) and U(x) so that (19) 
becomes applicable; for the one- loop case only the F(x). That is why we discuss here only 
the F(x). From (15), the F(x) may be written as a sum of N F < \N(N+1) non- vanishing, 
bilinear terms in xf. 



F( x yVf»-4L/2) = 

I ^ N F Np 

= m - dm) (2^f n_ ( y +w dz - n iMv^r 

lMij)XiXj] -^-^-^ r ( k - d ± + £ z) n r(-*j). 

(20) 

Here, f n (hj) = /ij if /ij 7^ 0. Inserting (20) (and if needed a similar representation for 
the U(x)) and the tensor function P(T) into (2) allows to apply (19) for an evaluation of 
the ^-integrations. 

As a result, any scalar Feynman integral may be represented by a single multi-dimensional 
MB-integral and L-loop tensor integrals by finite sums of MB-integrals. With AMBRE we 
will evaluate the L-loop integrals by a loop-by-loop technique, which essentially allows us 
to restrict the formalism to the one-loop case. By the examples it will be seen that this 
is a powerful ansatz for many applications. 

In subsequent steps, the package MB may be called. This package needs as input some 
MB-integral(s), e.g. as being prepared by AMBRE. As described in detail in [8], MB allows 
to analytically expand a Feynman integral in e and to evaluate the resulting sequence of 
finite MB-integrals by one or the other method. 



N F 



fn{h 



3 Using AMBRE 



In this section we describe the use of the Mathematica package AMBRE. AMBRE stands for 
Automatic Mellin-Barnes Representation. It is a (semi-)automatic procedure written for 
multi-loop calculations. The package works with Mathematica 5.0 and later versions of it. 

The algorithm to build up MB-representations for Feynman integrals as described in the 
last section consists of the following parts: 

(i) define kinematical invariants which depend on the external momenta; 
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(ii) make a decision about the order in which L 1-loop subloops (L > 1) will be worked 
out sequentially; 

(iii) construct a Feynman integral for the chosen subloop and perform manipulations on 
the corresponding F-polynomial to make it optimal for later use of the MB repre- 
sentations; 

(iv) use equation (20); 

(v) perform the integrations over Feynman parameters with equation (19); 

(vi) go back to step (iii) and repeat the steps for the next subloop until F in the last, L th 
subloop will be changed into an MB-integral. 



The steps (ii) and (iii) must be analyzed carefully, because there exists some freedom of 
choice on the order of loop integrations in step (ii) and also on the order of MB integrations 
in step (iii). Different choices may lead to different forms of MB- representations. 

The present version 1.0 of AMBRE can be used to construct planar Mellin-Barnes represen- 
tations for: 

• scalar multi-loop, multi-leg integrals 

• tensor one-loop integrals 

• integrals with specific higher-rank numerators ending up with a single MB-integral 

In the next sections several examples will be used for an introduction to specific features 
of the package. 

Here, we describe basic functions of the package. The starting point of all calculations is a 
proper definition of the integral (2) and of the kinematical invariants to be used. Formally, 
it has to be done in the following way: 

Fullintegral [{numerator} , -[denominator}, {internal momenta}]; 
invariants = {invariants as a rule}; (22) 

We recommend to use hi and Pi as symbols for internal and external momenta, respectively. 
Also non-zero masses should appear as symbols; a numeric value may cause problems in 
multi-loop calculations. 

The command Fullintegral defines a given integral. For example: 



Fullintegral [{1}, {PR[kl, 0, nl]*PR[kl + pi, m, n2]}, {kl}] ; (23) 



(21) 



corresponds to: 




(24) 
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The last argument in the Fullintegral function is a list of internal momenta. The order 
of internal momenta in this list controls the ordering of integrations (if iterated). For 
example {k3,k2,kl} defines the first integration to be over k 3 , the second over k 2 and 
the third over k x . The next step is to prepare a subloop of the full integral by collecting 
all propagators which carry a given loop momentum fcj. We do this by initiating the 
consecutive functions: 

IntPart [iteration] (25) 

Each iteration, % — 1, 2, . . . , L, prepares the appropriate subloop for the integration over 
the corresponding internal momentum. It will display a piece of the Fullintegral with: 

• the numerator associated with the given subloop 

• subloop for a given internal momentum 

• internal momentum for which AMBRE will integrate the subloop 

The execution of IntPart [iteration] proceeds in the order IntPart [1], IntPart [2] , 
then IntPart [3] , and so on. If there is a need to change the ordering of integrations, 
one has to change the order in the starting list of internal momenta (22). Inserting 
IntPart [2] before IntPart [1] would not be a proper way to do this. In the output 
of IntPart [iteration] a tag message will be displayed: 

Fauto : :mode : U and F polynomials will be calculated 

in AUTO mode . In order to use MANUAL mode execute Fauto [0] . (26) 

By running Fauto [0] , AMBRE will calculate the F-polynomial (with name f upc) for a 
given subloop. At this stage, a user might wish to modify f upc manually, e.g. by applying 
some changes in kinematics. 

During the calculations, the FX function of AMBRE may appear in the F-polynomial. This 
function collects full squares of sums of Feynman parameters, e.g.: 



Such terms appear in the F-polynomials if some masses in the loops are equal. They will 
later allow to apply Barnes' lemma leading to lower dimensional MB-representations. At 
the other hand, the exponent two of the square may lead to arguments of T-functions in 
(20) with doubled integration variables, with far-reaching consequences for an analytical 
evaluation when a sum over an infinite series of residua is tried. 

The basic function for deriving the Mellin-Barnes representation is: 

SubLoop [integral] (28) 




(27) 
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This function takes output generated by IntPart [iteration] and performs the following 
calculations: 

• calculate the F-polynomial for the subloop (only if Fauto [1] is set) 

• determine the MB-representation for the F-polynomial 

• integrate over Feynman parameters X{ 

As a result, the MB-representation for a given subloop integral will be displayed. In multi- 
loop calculations one will notice additional propagators (marked in red in the output of 
AMBRE) which appear from the intermediate F-polynomial (see section 5.1 for an instruc- 
tive example). 

As mentioned, AMBRE can construct Mellin-Barnes representations for general one loop 
tensor integrals. The procedure of calculating such cases is basically the same, with few 
minor differences. First of all, the numerator input must be defined. A one-loop box 
diagram with numerators (}i\'P\){k\'P2){k\'Pz) might look like this: 

Fullintegral [ 

{kl*pl,kl*p2,kl*p3}, 

{PR[kl,m,nl]PR[kl+pl ) ) n2]PR[kl+pl+p2 ) m,n3]PR[kl+p3 ) 0,n4]},{kl}] ; 

(29) 

We have written this procedure such that numerators consist of scalar products of internal 
and external momenta. In the calculations with tensors, the definitions of momentum 
flows in the subloops play a crucial role for the results and have to be controlled carefully. 
Another difference to scalar cases is the way how AMBRE displays results. Because they 
can be long, we decided to use a short notation. For example: 

{ARint [1] , ARint [2] , ARint [3] } (30) 

The result of the evaluation has to be understood as the sum of the elements, 
ARint [1] +ARint [2] +ARint [3] , 

where each ARint [i] is one of the resulting MB-integrals. By executing 
ARint [result , i] 

one may display the appropriate ARint [i] . The procedure uses the short notation by 
default, but it is also possible to use the option Result->True in order to force SubLoop 
to display the full result: 

SubLoop [integral ,Result->True] ; (31) 
Finally, we have also implemented Barnes' first lemma: 
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f , „, Vn/ , x r(a + c)r(a + d)r(& + c)r(6 + d) 

— ioo 

and Barnes' second lemma: 

J 3 r(a + z)T(b + z)T{c + g)T(d - g)T(e - z) 
J Z T(a + b + c + d + e + z) 

—ioo 

r(g + d)T{a + e)r(6 + d)T{b + e)r(c + d)T{c + e) 3 
r(a + b + d + e)r(a + c + d + e)r(6 + c + d + e) ' ^ 

The usage of Barnes' lemmas is simple; one has to execute: 

BarnesLemma [representation, i] (34) 



where i is 1 or 2 for the first or second Barnes' lemma, respectively. This function tries 
to apply the lemma on all integration variables Z{ of the MB-representation which do not 
appear in the exponents of kinematical invariants. It also searches in the exponents of 
kinematical invariants for pairs of integration variables. For example, be in one exponent 
the combination (z,\ + z?) and in another one the combination {—z\ — z-i). This might 
appear as dictated by the structure of equation (1). The automatic change z 1 — > z 1 — z 2 
eliminates z 2 in these exponents so that Barnes' lemma can be applied for z 2 . A comment 
will be displayed if the lemma was successfully applied. Barnes' first lemma is quite often 
applicable, while Barnes' second lemma applies only sporadically (see example 14. nb). 
The automatic change of variables may be switched off by calling shift [0] . 

In the appendix we list the Mathematica functions of AMBRE. 



4 One-loop integrals 



We will give a couple of examples starting with construction of MB-representations for 
the 1-loop Feynman integrals which are an important ingredient of the algorithm (21). 
Most of the cases considered in subsequent sections are connected with massless gauge 
theories or massive QED. 



4-1 Example: the pentagon diagram of massive QED 



Let us consider the one- loop five-point function shown in figure 1. The external momenta 
fulfill p\ = 0, pf = m 2 for the other particles, and the = (pi + pj) 2 are kinematical 
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Fig. 1. Massive QED pentagon diagram 

invariants of the process. If we naively use the FUPolynomial function of the MB package, 
we will get: 

U = xi + x 2 + x 3 + x 4 + x 5 , (35) 
F = m 2 x\ + 2m 2 XiX 3 — S15X1X3 + m 2 xl + 2m 2 XiX 4 — S23X1X4 + m 2 x 2 X4 

— S45X2X4 + 2m 2 x^Xi + m 2 x\ — 512X2X5 + rn 2 x 3 x^ — S34X3X5. (36) 

A simple counting of terms in the F-polynomial would prove that this leads to a twelve- 
dimensional MB-integral. Of course the terms in F can be grouped from the beginning 
and we will see in a minute that a five-fold MB-integral may be obtained; see also the 
sample file examplel .n 

First, propagators and kinematical invariants are defined: 

Fullintegral [ 
{1}, 

{PR[kl + pi, 1, nl]*PR[kl + pi + p5, 0, n2]* 

PR[kl + pi + p4 + p5, 1, n3]*PR[kl + pi + p3 + p4 + p5, 1, n4] * 
PR[kl + pi + p2 + p3 + p4 + p5, 0, n5]},{kl}]; (37) 

The kinematics is defined in a cyclic way: 

p 2 = m- , s i>i+1 = (pi + Pi+i) 2 , i = 1, 5. (38) 

Then, using the IntPart and SubLoop functions the steps (ii)-(v) of the algorithm (21) 
are worked out and we end up with a nine-fold MB-representation. This representation is 
due to the following F-polynomial, constructed in the automatic way by AMBRE: 

F = m 2 (x 1 + x 3 + x 4 ) 2 - S15X1X3 - S23X1X4 + m 2 x 2 X4 - s^x 2 x^ - Syi^h 

+ m 2 x 3 x 5 - S34X3X5. (39) 

2 The sample Mathematica files are part of the package AMBRE. They are available at [9,10]. 
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Some mass terms have been collected here, but the F-polynomial can be further simplified 
by redefining s 34 — > s 34 + m 2 and s 45 — > s 45 + m 2 , so that each term XiXj appears only 
once. The F polynomial becomes finally: 



F" = m 2 (xi + x 3 + X4) 2 - si 5 xix 3 - s 23 x 1 x i - s 45 x 2 x 4 - si 2 x 2 x 5 - s 34 x 3 a; 5 , (40) 

which gives a seven-fold MB-representation. In certain cases, some of the MB-integrations 
do not depend on the kinematics and Barnes lemmas may be applied. Here, due to the term 
(xi+a; 3 +:r 4 ) 2 one may twice use Barnes' first lemma (32) and thus the MB-representation 
can be further reduced to a five-fold integral. A five-particle scattering process depends on 
five variables (plus a mass in Bhabha scattering), so a further simplification is impossible. 

In sample file example2 .nb, we use another definition of kinematical variables, namely 



pi 



2 

m , 



P1P2 



P1P4 = m 2 + \(vi - s - t), pip 5 
P2P4 = \{s-vi-v 2 - 2m 2 ), p 2 p 5 



PzPa = 2 v 2 



|(t'-2m 2 ), p 4 p 3 
|(s-2m 2 ), p 2 p 3 
\{v 2 - s -f + 2m 2 ), 

PiPb 



(41) 



P3P5 = \{t' ~t-V 2 ), 



Ht-2m 2 ), 



and we get F directly in the form: 



F = (xi + x 3 + X4) 2 — txix% — t'x\X 4 — v 2 x 2 x 4 — sx 2 x 5 — Vix 3 x 5 . (42) 

No wonder, that using function SubLoop we obtain directly the smallest, seven-dimensional 
integral, which then again reduces to the five-fold integral. The resulting MB-representa- 
tion for the scalar Feynman integral is: 



g -i-ioo+Ui 

°n = w n / *.(-«)-- "nr(-o- (^)" 4 (^y 

l - L -ioo+Ui 

r[-r 2 ]r[-r 3 ]r[i + r 2 + r 3 ]r[-n + r 2 + r 3 ]r[-2 - e - n - r 4 ]r[-r 4 ] 
r[l + r 2 + r 4 ]r[-2 - e - n - r 5 ]r[-r 5 ]r[l + r 3 + r 5 ] 

r[3 + £ + n + r 4 + r 5 ]r[3 + 2ri + r 4 + rj f , 

r[-l-2e]r[3 + 2(r 2 + r 3 )+r 4 + r 5 ] ' 1 j 

The real parts of the integration strips are — 2 < ui < — 1 and — 1/2 < < 0, i = 2 . . . 5. 

A subsequent application of MB shows that up to constant terms in e, needed for an evalu- 
ation of two-loop massive Bhabha scattering [11], there are maximally three-dimensional 
finite contributions to be evaluated further. 



13 



Fig. 2. Massive QED one-loop box diagram 

4-2 Numerators 

AMBRE may handle arbitrary one-loop tensor integrals. The one-loop Feynman parameter 
integral for a tensor of degree m is the generalization of equation (18): 

G 1 (T m )=G 1 (k^---k^) 

(44) 

Here F = F(x) and P = Px(A; M ) = J2i x iPi — J2i.e x iPiePe were introduced in equations 
(9) and (17). The r starts from zero (with Aq = 1), and it is A r = for r odd, and 
•Ar = g^ 1 ^ 2 ■ ■ ■ g^r-i^A f or r even. The convention . . .] means the totally symmetric 
combination of the arguments. 

In AMBRE tensorial numerators are assumed to be contracted with the external momenta 
p e , so that the following quantity is evaluated: 

P m G 1 (T m ) = (p^...p^) G 1 (k^---k^). (45) 

As an example, we have prepared the massive QED one-loop box of figure 2 in sample file 
example3.nb with the numerator {kiPx){k\P2){k\Pz). The corresponding definition used 
in AMBRE is: 

Fullintegral [ 

{kl*pl,kl*p2,kl*p3}, 

{PR [kl , m , nl] PR [kl+pl , , n2] PR [kl+pl+p2 , m , n3] PR [kl+p3 , , n4] } , {kl}] ; 

(46) 

Obviously, when working with tensor integrals we expect the result to be a sum of several 
MB-integrals (the higher the rank is, the more integrals will be obtained). We have cross 
checked numerically results for two-, three- and four-point functions by comparing our 
results (from using AMBRE and MB) with decompositions of integrals into master integrals 
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Fig. 3. General one- loop vertex 

using the Integration-By-Parts method implemented in the package IdSolver (M. Cza- 
kon, unpublished). Cross checks were done for numerators with up to eight scalar products 
in the numerators of the Feynman integrals. 

Finally we refer to section 5.2 for the interesting special case of irreducible numerators 
arising in intermediate subloops. In certain cases, the result for a tensor integral may 
remain as compact as it is for scalar integrals. 



4-3 More masses 



iV-point functions with arbitrary internal masses and off-shell external legs give compli- 
cated multi-dimensional MB-integrals. Let us consider here and in example4.nb a general 
one-loop scalar vertex, Fig. 3. In this case we get a five-dimensional MB-integral: 



(-If 123 

^general 



r\-zA 



T[n 1 ]T[n 2 ]T[n 3 ]T[A-2e-n 123 ] 

^ +200 +ioo +200 +ioo +ioo ^ 

— / dz\ \ dz 2 / dz 3 / dzd / dz 5 TT — — 

(2vri) 5 J J J A J J 5 f}, T 2 - e - n x - z x + z 5 ] 

—ioo —ioo —ioo —ioo —ioo 

(m^ (^f (™D 2 ~™™ 34 5 (ml + mj- Ml) z \m\ + ml - M 2 2 ) 24 

(m 2 2 + m 2 3 - Ml - M 2 2 - s) Z5 Y[ ni + z 1124 ]r[4 - 2e - n U23 - z U24 \ 

r[n 2 + z 235 ]T[-2 + e + n l23 + ^2345], (47) 

where we abbreviated ^124 = 2z\ + z 2 + z 4 and n 123 = rii + n 2 + n 3 , etc. 

For the massive QED case, M\ = M 2 = m 2 = m 3 = m, mi = 0, we get a compact 
one-dimensional MB-representation: 



_l)»i»r[4-2e- mi23 ] 1 , ., 2 .,, .,-,-„ m -. 



^ ED ^ r[, i2 ]rMr[4-2e- , il23 ] ^ / dz{m > { ~ s 

—ioo 

T[-z}r[2 - e - ma ~ z]T[2 - e - n 13 - z]T[-2 + e + n 123 + z] 

r[4 - 2e - m 123 - 2z] ' 1 
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Fig. 4. Six-point scalar functions; left: massless case, right: massive case 



4-4 More legs 



For topologies with a higher number of legs, there is an increasing number of kinematical 
invariants and so the dimension of MB-representations increases. The number of dimen- 
sions may become smaller after analytical continuation in e and for lowest orders in e. 
For a scalar or vector Bhabha massive five-point function, Fig. 1, up to constant terms 
in e, it includes at most three-dimensional integrals, which hopefully can be solved even 
analytically [11]. In general, the MB- representation for that case is five-dimensional, see 
section 1. 

In example5.nb we derive MB-representations for a massless and a massive one-loop 
hexagon scalar diagram, see figure 4. In general, it is an eight-fold integral, but the constant 
term in e includes again only up to three-dimensional MB-integrals. 

If all internal lines have equal non- vanishing mass, one has to deal with a nine-dimensional 
MB-integral. Again, the numerical results have been checked for both cases in the Eu- 
clidean region against sector decomposition. The package contains the auxiliary file 
KinematicsGen.m which generates the kinematics for six-point functions with arbitrary 
off-shell external legs. 



5 Multi-loop integrals: loop-by-loop integrations 



The Feynman integral (5) includes a delta-function which makes U = 1 for one-loop 
diagrams there so that the MB-relation (1) acts only on F. This simplification can be made 
also useful in multi-loop integrals by performing loop-by-loop integrations. We collected 
few examples which will exhibit several specific features. 
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Fig. 5. Massive two-loop planar QED box 
5. 1 Example: two-loop planar box in massive QED 

Let us take the massive two-loop planar box topologjOO with seven internal lines as in- 
troduced in example6 .nb. The momentum flow is defined in the following way, with all 
momenta being incoming: 

Fullintegral [ 
{1}, 

{PR[kl, m, nl]PR[kl + pi, 0, n2]PR[kl + pi + p2, m, n3] 

PR[kl - k2, 0, n4]PR[k2, m, n5]PR[k2 + pi + p2, m, n.6] 

PR[k2 - p3, 0, n7]}, {k2, kl}] . (49) 

First, the momentum integration over hi is taken. The ki flow in the first subloop is 
defined by the function IntPart [1] , which contains all propagators with momentum k%: 

integral = PR[kl - k2, 0, n4]*PR[k2, m, n5] * 

PR[k2 + pi + p2, m,n6]*PR[k2 - p3, 0, n7] . (50) 

We just mention that generally it is preferred to choose the order of iteration such that 
first the loops with lowest number of lines are executed. Then their F-polynomials have 
a minimal number of terms. The first loop's F-polynomial is the SubLoop [integral] 
function: 

F[k 2 ] = fupc= 

m"2*FX[X[2] + X[3]]~2 - PR[kl, m] *X [1] *X [2] - 

PR[kl + pi + p2, m]*X[l]*X[3] - s*X[2]*X[3] - 

PR[kl - p3, 0]*X[1]*X[4] + 4*nT2*X[3]*X[4] - 

s*X[3]*X[4] - t*X[3]*X[4] - u*X[3]*X[4] (51) 

It is reproduced here as derived without interactions by the user. The F-polnomial con- 
tains a mass term with the FX function which later will allow to apply Barnes' first 

3 In fact there are three double-box diagrams in massive QED. One of them is non-planar, and 
we discuss here the so-called first planar diagram [3]. 



17 



lemma successfully, and also a redundancy in X [3] *X [4] . The following nine-fold Mis- 
representation after integrating over k 2 is obtained: 

SubLoopl[((-l))~(n4 + n5 + n6 + n7 + z2 + z3 + z5) 4~z6 
(m~2~(zl + z6) (-s)-(z4 + z7) (-t)~z8 

(-u)~(2-ep - n4-n5-n6-n7 - zl - z2 - z3 - z4 - z5 - z6 - z7 - z8) 

Gamma [-zl] Gamma [(-z2)] Gamma [(-z3)] Gamma [ 

2 - ep - n4 - n5 - n6 - zl - z2 - z3 - z4] Gamma [(-z4)] 

Gamma [(-z5)] Gamma [n4 + z2 + z3 + z5] Gamma [(-z6)] Gamma [(-z7)] 

Gamma [-z8] Gamma [-2+ep + n4 + n5 + n6 + n7 + zl + z2 + z3 + z4 + z5 + 

z6 + z7 + z8] Gamma [ 2 - ep - n4 - n5 - n7 + zl - z2 - z5 - z9] 

Gamma[(-z9)] Gamma[(-2) zl + z9] Gamma[n5 + z2 + z4 + z9])/ 

(Gamma [n4] Gamma [n5] Gamma [n6] Gamma [4 - 2 ep - n4 - n5 - n6 - n7] 

Gamma [n7] Gamma [(-2) zl])), 

PR[kl, m, z2]PR[kl + pi + p2, m, z3]PR[kl - p3, 0, z5]]) (52) 

It is clear that the factors in front of the X[3]X[4] coefficient sum up to zero, due to 
s + t + u = 4m 2 . To remove them from the beginning, the Fauto [0] option must be 
executed, followed by a modification of F: 

fupc = fupc /. u -> 4*m~2-s-t. (53) 

In this way, executing the SubLoop [integral] function again, the MB-representation 
becomes five-dimensional, and also the term 4 Z6 is absent now. 

The same situation appears in the second iteration, when integrating over k\. We can 
switch to the Fauto [0] mode and again modify F. After again applying Barnes' first 
lemma, we end up with a six-dimensional integral. 

Of course, by writing from the very beginning the kinematical invariants without the 
invariant u, one can work out the whole case fully automatic with mode Fauto [1] . 

5.2 Special numerators 

The example6.nb is interesting in yet another respect. After the first integration, the 
propagators for the second one contain four propagators, some of them with shifted indices 
compared to the input: 

PR[kl, m, nl-z2]PR[kl + pi, 0, n2]PR[kl + pi + p2, m, n3-z3] 

PR[kl - p3, 0, z5]. (54) 
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This corresponds to the one-loop box of example example3 .nb discussed in section 4.2, but 
with shifted indices. It includes the one new propagator with momentum q 5 = k\ —p 3 . If we 
would have been evaluating an integral with numerator (gf )~ ng and repeat the calculation, 
we would get after the k 2 integral an F-polynomial with one of the terms including the 
propagator PR[kl + pi + p2 + p4, 0, 1] ; see SubLoop [integral] in example7.nb; 
see also [5]. It will sum up with PR[kl + pi + p2 + p4, 0, -n8] resulting in the fol- 
lowing integral 

integral= PR[kl, m, nl - z2] PR[kl + pi, 0, n2] 
PR[kl + pi + p2, m, n3 - z3] 
PR[kl + pi + p2 + p4, 0, -2 + ep + n4 + 

n5 + n6 + n7 - n8 + zl + z2 + z3 + z4] , (55) 

which has the following well-known F-form of the one-loop box: 

m~2 FX[X[1] + X[3]]"2 - s X[l] X[3] - t X[2] X[4]. 

What is essential here, no additional momentum structure appears. 

Analyzing the irreducible numerators of the topology for the given momentum choice, 
one finds that there are two scalar products which may not be represented by linear 
combinations of the propagators (and thus are called irreducible): kip 3 together with k 2 p± 
or kip 3 together with k 2 p 2 . So, q\ represents one of two existing irreducible numerators 
and it may be quite useful to have a simple MB-representation for that case. We see 
that there are integrals with (selected) numerators which may be represented by a single 
MB-representation as if a scalar integral would have been studied. This was used several 
times in examples given in [5,6] and in references cited therein, and it was also used e.g. 
in [12] for a study of massive two-loop box master integrals, and for more sophisticated 
four- loop cases in [13]. 

Finally, a six-dimensional MB-integral emerges like in the scalar case. To check this inte- 
gral numerically with the MB package, two analytical continuations, one in e and one in 
one of the powers of propagators must be done. We have checked the numerical result also 
against the results we got from a sector decomposition calculation and from a small-mass 
expanded version [14,12]. 



5.3 Further examples: A three-loop planar box, a four-loop self-energy, and a two-loop 
pentagon 

A three-loop planar integral, shown in figure 6, is treated in example8 .nb. The result is 
a 10-fold MB-representation. With the MB package it was shown that the numerical result 
agrees with [6] . 
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Fig. 6. The loop-by-loop iterative procedure 



The dimensions of some MB-representations for several massless and massive ladder 
topologies are summarized in Table 1. We apply an iterative procedure. For planar topolo- 
gies the loop-by-loop iteration gives always proper topologies which obey momentum con- 
servation. Only some powers of propagators change into non-integer (complex) numbers. 

A similar procedure can be applied to more complicated topologies which obey the same 
rule: integrating over an internal momentum leads to a topology with propagators and mo- 
mentum flow obeying momentum conservation in the remaining parts, i.e. we get regular 
subtopologies. 

In this procedure, the choice of momenta flowing and the order of iterations are very im- 
portant. Look e.g. at the two-loop ladder example, also shown in figure 6. If we would allow 
for the momentum flow k\ through all the outer lines, and take first the integration over k\ 
and then that over k 2 , the final representation would not come out optimal (and Barnes' 
lemmas do not help). Starting instead with the k 2 integration, we will again end up, as 
with the momentum flows shown in the figure, with a six- dimensional representation. 

In files example9 .nb and examplelO .nb, massless MB-representations are constructed for 
a four-loop two-point topology and for a two-loop five-point massless topology, see figure 7. 
The six- dimensional four-loop self-energy has been checked numerically against sector de- 
composition. In examplelO .nb, there are three different derivations of MB-representations 
for the same kinematics, defined by equation (38). In each case we got another dimension 
of MB-integrals. The minimal dimension of the integral is seven when we integrate first 
over internal momenta of the box and then over that of the pentagon. We checked that 
this agrees numerically with [15] where also a seven- dimensional MB-integral has been 
obtained. If we integrate first over the internal momentum running in the pentagon and 
next over that in the box, then a nine- dimensional MB-integral is obtained; again numeri- 
cally they agree. In the third derivation, the momentum flow in the propagators is chosen 



Massless 


Massive 


1-loop 


2-loop 


3-loop 


4-loop 


1-loop 


2-loop 


3-loop 


4-loop 


1 


4 


7 


10 


3 


8 


13 


18 


1 


4 


7 


10 


2 


6 


10 


14 



Table 1 



Dimensions of ladder topologies before and after applying Barnes' first Lemma. 
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Fig. 7. Massless topologies; left: four-loop two-point diagram, right: two-loop five-point diagram 
in a different way. Then a 13- dimensional MB-integral results. 



6 Tadpoles 



The loop-by-loop approach can also be applied to planar tadpoles. Attention must be paid 
to keep the right order of integrations. Making iterations with the Fauto [1] option (i.e. 
automatic), we may end up with three different forms of propagators in the last iteration: 
one massive propagator, massive and massless propagators, or one massless propagator. 
For the first situation the well known formula is used in AMBRE: 

' ^^(-rr^,- 2 ',,, 1 . .. (56) 



(k 2 -q 2 f v ' Y\v\ (q 2 )^ 



e-2 



We found that for some massive tadpoles a term (— m) a can appear which would lead to 
an oscillatory error while doing numerical calculations with MB. In such a situation one 
has to go back to the previous subloop and modify the F-polynomial with Fauto [0] so 
that two propagators with equal momenta appear: a massive and a massless one. The 
same procedure must be applied when a single massless propagator appears in the last 
integral. 

We give as an example examplell .nb, for the diagram also shown left in figure 8. Using 
AMBRE, we have constructed a one-dimensional Mellin-Barnes representation: 



r it i +200 

T(n u ...,».) = (m^-^ ^ 1 " ' -Vff f ' " / *» 

—ioo 

x r[-zi]r[2 - e - m - zi]r[2 - e - n 2 - Z!]r[4 - 2e - n 12 - z x \ 
r[-6 + 3e + ni245 + Zi]r[-8 + 4e + ni 234 5 + z x ) 



T[4 - 2e - ma - 2«i] 



(57) 



At this point the proper order of integrations is very important. A different choice can 
lead to two- or even higher-dimensional representations. 

Using MB we got for the basic integral numerically: 
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Fig. 8. Four-loop tadpoles with three massive lines; left: with a one-dimensional MB-represen- 
tation, right: with a six-dimensional MB-representation 

T(l, 1, 1, 1, 1) = -4^ + 4 + 2.8433003667574474 + 5.781543610421033- 

e 4 e d e 2 e 

+ 22.955621881705923 + 80.89550616785341e + 1085.2836587072804e 2 
+ 4545.303884134432e 3 + 35998.99383263255e 4 , (58) 

This is in agreement with [16]. 

However, it appears that MB-representations for four-loop tadpoles can be more compli- 
cated. In examplel2 .nb, treating the diagram in figure 8 (right), we get a six- dimensional 
MB-integral. Taking into account other approaches [16,17], one may see that the MB- 
approach to multi-loop calculations has natural limits, especially in the massive cases. 



7 On-shell diagrams 



Mellin-Barnes representations can be also useful for solving on-shell topologies. For on- 
shell self-energies, one may use the package 0N-SHELL2 [18,19] written in FORM v.2.3 
[20]. In examplel3.nb we show how to evaluate the self-energy SE513m2 shown in figure 
9, which is in the notations of [18] the diagram F01101. The MB-representation is a 
two-dimensional integral: 



(27ri) 2 r[l -2e] J J V[l - 3e - Zl \ 



—too — iOO 



x r[-e - z_i - z 2 ]r[-3e -zi + z 2 ]T[l - e + z x + z 2 ]T[l + e + z x + z 2 ] , . 

T[l-e-z 1 + z 2 ]T[2 + e + z 1 + z 2 ] ' 1 ' 

In the example, the agreement with the result of 0n-Shell2 result is demonstrated. 

A simpler case is SE311m with one massive and two massless propagators, see figure 9. 
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e 



Fig. 9. On-shell self-energies; left: SE513m2, right: SE311m 



The result is simple: 



SE311mOS = -(m 2 ) (1 ~ 2e) 



r[3-4e]r[l -e] 2 r[e]r[-l + 2e] 
r[3-3e]r[2-2e] 



(60) 



It can be expanded easily to any order in e. Here, in example 14 . nb, Barnes' second lemma 
has been used, which happens not too often. Again, the agreement with the 0n-Shell2 
result is presented. 



8 Non-planar topologies 

The loop-by-loop iterative procedure described in this paper seems to be not the most 
efficient approach in the case of non-planar topologies. It is known from [21] that the 
massless non-planar vertex is described by a two-dimensional Feynman parameter integral. 
If we consider the loop-by-loop procedure for this case, we can divide the two-loop topology 
in figure 10 into two parts (follow the vertical line). The hourglass topology on the right- 
hand side, with two off-shell legs, gives a three-dimensional MB-representation [6], and 
adding the second part on the left-hand side we end up with a four-dimensional integral. 
No matter how we arrange the momenta flows in the diagram, it cannot become better. 
To get the minimal, two-dimensional integral, another approach must be realized. It is an 
open question to us if the representation of non-planar diagrams can be automatized in 
a way like that for planar caseJ^l 



The non-planar examples in the study [7] do not go beyond our observations stated here. 




Fig. 10. Non-planar massless vertex 
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9 Summary 



We have described the Mathematica package AMBRE for the construction of MB-represen- 
tations for planar Feynman integrals and gave in a tutorial part a variety of sample 
applications. Typically, the iterative loop-by-loop approach gives a possibility to construct 
MB-integrals of minimal dimension. Usually Barnes' first and second lemmas help to 
get the minimal dimension of MB-integrals, independent both of the flow of momenta 
in diagrams and of the order of iterations. However, for more complicated kinematics, 
starting with five legs, the order of iterations and the choice of momenta flows matters. 
As is shown there in the case of tadpoles, MB-representations for massive topologies are 
not always the best way of evaluation. For some topologies quite simple representations 
are found, however, also multi-dimensional MB-integrals may arise from which it is hard 
to get stable, accurate numerical results, not mentioning exact analytical results. 

Constructing useful MB-representations for a given Feynman integral is a kind of an art. 
As an example, let us mention the QED master integral B512m2 (a diagram with five lines, 
two of them being massive; notations are due to [22,12]). This integral may be obtained 
by contracting directly two lines in the massive Bhabha two-loop planar integral B714ml 
[3] (the so-called first planar master of massive QED). In [12] it was shown that, after 
expansion in e, the expression for B512m2 consists of eleven integrals, one being four- 
dimensional. This was compared to constructing B512m2 from the scratch, loop-by-loop. 
Here, again after expansion in e, we are left with four integrals, all of them being three- 
dimensional or simpler. This can be checked easily by the reader using the MB package and 
both the representation B714ml given in [3] (the contraction of two lines must be done 
there) and the representation B512m2 given in [12]. There is no simple relation between 
both representations, due to our lack of knowledge on more complicated relations between 
integrals of different dimensionality. One should also mention that an independent check of 
the MB-representations always is strongly recommended. By construction, neither AMBRE 
nor MB perform rigorous proofs of their applicability. 

Certainly, the number of integration variables is of importance for the final evaluation of 
MB-integrals, both in a fully analytical form or using approximations in some kinematical 
limits. In many cases some package like XSUMMER [23,24] can be used after deriving sums 
over residua. This again might become non-effective if the number of the nested sums - 
connected with the dimension of the MB-integrals - is too large or if the result is not in 
the class of functions covered by (e.g.) XSUMMER. Similar statements hold for the case of a 
fully numerical evaluation of MB-integrals. 

To summarize, for many applications of present phenomenological or more theoretical 
interest the package AMBRE solves an important part of the complete calculational problem: 
the derivation of expressions for a large class of Feynman integrals, which may then be 
used for further study. 
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A AMBRE functions list 

The basic functions of AMBRE are: 

• Fullintegral[{numerator}, {propagators}, {internal momenta}] is the basic 
function for input Feynman integrals 

• invariants is a list of invariants, e.g. invariants = {pl*pl — > s} 

• IntPart [iteration] - prepares a subintegral for a given internal momentum by collect- 
ing the related numerator, propagators, integration momentum 

• Subloop [integral] - determines for the selected subintegral the U and F polynomials 
and an MB-representation 

• ARint [result, i_] - displays the MB-representation number i for Feynman integrals 
with numerators 

• Fauto[0] - allows user specified modifications of the F polynomial fupc 

• BarnesLemma[repr,l,Shifts->True] - function tries to apply Barnes' first lemma 
to a given MB-representation; when Shifts->True is set, AMBRE will try a simplifying 
shift of variables 

BarnesLemma[repr,2,Shifts->True] - function tries to apply Barnes' second lemma 
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